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Abstract 

^h ; 

Exclusive vector meson production reactions such as Au + Au — » Au + 

Mh ■ Au + V, where V = p,ui,<p or J/ip can proceed through photon- Pomeron 

\& , and photon-meson interactions. Photons from the electromagnetic field of 

one nucleus interact coherently with the other nucleus. Photonuclear cross 

^- ' sections are scaled from jp data, and convoluted with the photon spectrum 

0\ ■ to find the exclusive rates. 

^vq ■ The cross sections at the RHIC and LHC heavy ion colliders are huge, 

10% of the total hadronic cross section at RHIC, and 50% at LHC. These 
qv | accelerators may be useful as vector meson factories. With iodine beams at 

ON ' RHIC, 640 p are produced each second (10 10 /year); with calcium at the LHC 

the rate is 240 kHz. The cj) rates are 39 Hz at RHIC and 15 kHz at LHC, 
while the J/ip rate is 0.3 Hz at RHIC and 780 Hz at the LHC. 
Qh. Because of the coherent couplings, the reactions kinematics are similar to 

r~| ! coherent two-photon interactions; we discuss the interplay between the two 

^| | reactions. 

J^J ■ (Submitted to Physical Review C.) 
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I. INTRODUCTION 

The electromagnetic fields of ultrarelativistic heavy ions may be treated as fields of 
virtual photons. When two nuclei pass by at small impact parameters, the photon field 
of one nucleus can produce a photonuclear interaction in the other nucleus. Most studies 
of photonuclear reactions in heavy ion collisions have considered hard interactions, such as 
heavy quark production |lj or photonuclear breakup p. 

We consider here the exclusive reaction A + A — > A + A + V, where A is a heavy nucleus, 
and V is a vector meson. This reaction can proceed via photon-Pomeron or photon-meson 
interactions, where the photon comes from the electromagnetic field of one nucleus, and the 
Pomeron or meson couples to the other nucleus. At sufficiently high energies, the momentum 
transfers from the nuclei are small enough that both couplings are coherent to the entire 
nucleus. Then, the nucleus remains intact, or nearly intact. We calculate the production 
rates and rapidity distributions for the p°, u, 4> and J/-0. 

These reactions will be studied experimentally at the upcoming Relativistic Heavy Ion 
Collider (RHIC) || now under construction at Brookhaven National Laboratory [|J, and 
the Large Hadron Collider (LHC) |J being built at CERN. For later reference, Table | gives 
the expected luminosities for various beam species at these two colliders. 

II. PHOTON FLUX AND PHOTONUCLEAR CROSS SECTIONS 

The calculations are simplest in the rest frame of one of the nuclei, labelled the tar- 
get nucleus. The total production is found by taking each nucleus as the target in turn. 
This calculation divides into two parts: determination of the photonuclear interaction cross 
section, and calculation of the effective photon flux from the emitting nucleus. 

In this frame, the photon energy k is large, and the momentum transfer squared from 
the target nucleus is small. The momentum transfer t has longitudinal tu and transverse 
components t±, with t = t\\ +t±. Since the initial state longitudinal momentum is small, t± 
is the final state momentum transverse to the photon direction. When t± — 0, \t\ — \t m i n \ = 



M 2 /4k. 

The coherent coupling of the photons and Pomerons limits their virtuality to 1/Ra, 
where Ra is the nuclear radius, or about 30 MeV/c for heavy nuclei like gold or lead. This 
virtuality is small enough that it can be neglected here, with the photon effectively real in 
the target frame. 

The photon flux from a relativistic heavy nucleus is given by the Weizsacker- Williams 
approach. Here 7 = Ea/tua ^> 1, where Ea and uia are the energy and mass of the nucleus. 
The photon flux at a distance r from the nucleus is 



d 3 JV 7 (fc,r) Z 2 ax 2 T ,> 



dkd 2 r TT 2 kr 2 



K((x) (1) 



where k is the photon energy, Z is the nuclear charge, K\ a modified Bessel function with 
x = kr/j. The target frame 7 corresponds to 7 = 2T 2 — 1, where T is the Lorentz boost 
in the center of mass frame (laboratory frame for collider geometry). Here, and throughout 
the paper, we take Ti = c = 1. 



In an exclusive interaction, there can be no accompanying hadronic interactions. In a 
hard sphere model, this occurs if the impact parameter b > 2R A , where R A = 1.2A 1 / 3 
fm. We use a more accurate approach, calculating, the probability of one or more hadronic 
interactions as a function of impact parameter. The nuclear density for a nucleus A a distance 
s from its center is modelled with a Woods- Saxon distribution for symmetric nuclei 

PA ^ S) = l + exp[(s-R ws )/d] (2) 

where the radius Rws and skin depth d are based on fits to electron scattering data || and 
p° is fixed by the normalization. This electron scattering data probes the proton distribution 
in the nucleus. If the neutron are differently distributed from the protons, this could affect 
the cross section.. 

The interaction probability at a given impact parameter is related to the overlap function 

T AA (\b\) = Jd 2 fT A (f)T A (f- b) (3) 

where r and b are 2-dimensional vectors perpendicular to the direction of travel z. The 
nuclear thickness function is 

T A (f) = jdz PA {ylW+#) ( 4 ) 

The number of nucleon-nucleon collisions follows a Poisson distribution with mean 
T AA (b)aNN, with o"at at the total nucleon-nucleon interaction cross section. The probabil- 
ity of having no hadronic interactions is then 

PoH{b) = exp(-T AA (b)a NN ). (5) 

We use (Inn = 52 mb at a beam energy of 100 GeV, rising to 88 mb at 2.9 TeV fTJ. 

We assume that the range of the Pomeron or meson field is much smaller than the size 
of the nucleus [f|. For b < 2R A , a form factor can be required to accurately model the 
electromagnetic fields inside the emitting nucleus; we neglect this tiny correction. 

The total photon flux is given by Eq. (|l|) integrated over all r, modulated by the 
noninteraction probability. 

dNJk) f°° n ,„„ .,, f R rdr f 2n , t d z NJk, b + r coa(d))) 

/ 2nbdbP 0H (b) I — T I d<P 7l ' W,) (6) 
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The b integral runs over impact parameter and r and over the surface of the target 
nucleus, for a given b. This process averages the field strength over the nucleus, neglecting 
the variation with r. Since the vector meson production is coherent over the entire nucleus, 
this averaging is appropriate. The integral is evaluated numerically; the result is shown in 
Fig. [l|. This flux is approximately equal to the photon flux in the region r > 2R A , which 
can be found analytically: 

dN,(k) _ 2Z* a ( _ £ _ v 



dk irk \ 2 



where X = IKj^Jkj^. This approximation is shown by the dotted line in Fig. |]. 

The 7 A — > V A cross sections are found with a Glauber calculation, using the observed 
7P -^ Vp cross sections as input. Data on 7p — > Vp covers a wide energy range, from low 
energy fixed target photon beams |J to ep collisions at HERA |10| . These cross sections 
may be parameterized as 

da(-yp — > Vp) 



dt 



= b v (XW e + YW- V ) (8) 

i=0 

where W is the center of mass energy, and Table [IJ gives the constants by, X, Y, e and rj, 
which are determined from fits to data [[LI]]. The X, e term represents the Pomeron exchange 
portion of the cross section. It rises slowly with energy (e ~ 0.22). This slow rise has been 
attributed to the increasing size and density of the proton; the increasing W probes smaller 



proton momentum fractions x. The Y, r\ term, is for meson exchange, primarily the /o |~1 



The meson exchange term falls rapidly as W rises. For the and J/tp, meson exchange is 
heavily suppressed, and the reaction occurs only through Pomeron exchange. For the J/ip, 
e ^> 0.22; the rapid rise has been ascribed to either a breakdown of the soft Pomeron model 
or to threshold effects ||10|| . Because of this behavior, the J/ip calculations must be treated 
with caution, especially at center of mass energies beyond that probed by HERA. 

To relate this to nuclei, we make use of the optical theorem and an Eikonalization 
technique. Following vector dominance |§, 



dai^p — > Vp) 



dt 



Ati a da{Vp — > Vp) 



t=o fl dt 



(9) 



with a similar relationship true for nuclei. Here, t is the squared 4-momentum transfer 
between the proton and vector meson, a is the electromagnetic coupling constant, e 2 /he, 
and f v is the vector meson-photon coupling, 

with M v the vector meson mass and T v ^ e + e - the leptonic decay partial width. However, 
measurements of f v from leptonic decay widths disagree with photoproduction data |TB 



A correction is required to account for non-diagonal coupling through higher mass vector 
mesons, as implemented in the generalized vector dominance model (GVDM). For the p, 
GVDM increases f% by 20% |D|. Because / 2 is used twice, the error largely cancels here. 
Values for / 2 /47r are given in Table [IT]. 

Using the optical theorem, the total cross section is 



2 nr ^ ,« da(Vp^Vp) 



afjVp) = 167T- 



dt 



t=o 



(11) 



At W = 10 GeV, we find 24 mb, 26 mb, 12 mb, and 1.0 mb for the p, u, <p, and J/ip total 
cross sections respectively. These values are consistent with theoretical and experimental 
expectations 0. The scattering cross section from heavy nuclei can be found by a Glauber 
calculation 

a tot (VA) = fd 2 r(l - e -— (^) t aa«Y (12) 



The optical theorem for nucleus A and GVDM are then used to find 



da(jA -> VA) 



dt 



a°UVA) 



« " -?P (13) 



These <i<r/(it|t=o are shown in Fig. |3] for gold nuclei. Also shown are the results of scaling 
the 7p cross sections by A 2 and A A / 3 ||. The p, to, and results vary only slowly with 
energy. The minimum around photon-proton center of mass energy W 1V = 10 GeV for the 
p and to marks the transition from meson dominated to Pomeron dominated reactions. <fi 
production, which is mediated only by Pomerons, shows a monotonic energy dependence. 

For the light mesons, da/dt\t=o is almost independent of energy, flatter than the cor- 
responding proton functions, which appears unchanged in the A 2 and A 4 ' 3 scaling curves. 
In contrast, da/dt\t=o for the J/ip rises rapidly with energy. Part of the difference is the 
different 7p cross section energy dependence, and part is the smaller J/ip interaction cross 
section. Because a P N is large, p interact near the nuclear surface, and a pA is the geomet- 
ric cross section ttR 2 a which is almost independent of k. However, ctj/^n is much smaller, 
and the J/ip interactions occur throughout the nucleus. As long as ctj/^nRaPo <^ 1, the 
rising (Jj/^n will reflect itself in aj/^ A . These cases correspond to the black disk and weak 
absorption limits respectively. 

The HERA results which provide the high energy 7p input data extend up to about W 1P 
= 180 GeV (k = 16 TeV in the target frame), with higher energy points extrapolations. 
Although the HERA data covers the entire RHIC energy range, the extrapolations are 
required for LHC. The uncertainty is especially problematic for the J/ip, where da/dt\ t= o 
must eventually flatten out, to avoid dominating the total photonuclear cross section. 

The total cross section depends on the slope of da /dt. For proton targets, the cross 
section can be parameterized as da /dt = A v exp(— bt + c\t\ 2 ) ||T^|. The slope b, and derivative 
c represent the sum of the proton size and the range of the interaction. Here, for simplicity 
we take c = with little loss of accuracy. Nuclei are much larger than protons, so b is 
dominated by the nuclear size, with b ~ R 2 and da /dt is dominated by the nuclear form 
factor. This is important because experimental study of these exclusive reactions depends 
on their small t± |15| which is determined by the form factor. 

The form factor is the Fourier transform of the nuclear density profile. Unfortunately, 
the Woods-Saxon distribution does not have an analytic form factor. Instead, for simplicity 
we approximate the Woods-Saxon distribution as a hard sphere, with radius R A , convo- 
luted with a Yukawa potential with range r = 0.7 fm [ IE|| . The Fourier transform of this 



convolution is the product of the two individual transforms: 

47rp, 



F(q=J\t\) 



Aq 3 



sm(qR a ) - qR a cos(qR A ) 



1 + a 2 q 



2„2 



(14) 



Figure |2] compares Eq. ([14]) with the numerical transform of the Woods-Saxon distribution 
for gold. The agreement is excellent. 

In addition to coherent production, there is also incoherent vector meson production, 
where the photon interacts with a single nucleon in the target. Incoherent interactions have 
a larger average t, and occur at a somewhat lower rate. They will not be further considered 
here, but they will be an essential component of any experimental analysis. 

The photonuclear cross section is 



a(jA -> VA) 



da (7 A -► VA) 



dt 



4=0 Jt 



dt\F{t)\- 



(15) 



For narrow resonances, t min = (M 2 /4k) 2 . Because of its width, the p is more complicated, 
and the cross section must be calculated using a Breit-Wigner resonance: 



da 
dMn 



a T p M p M n 



* (C-^F + rw' 



(16) 



Here, M p is the pole position of the resonance, a the total cross section (neglecting phase 
space corrections) and M n7T the final state invariant mass. The Breit-Wigner width T p 
includes the phase space correction, with T p = Tq^/pqYIMp/M^) [|17| where p n is the 
decay pion momentum in the p rest frame, with po= 358 MeV/c the pion momentum for 
M nn = M p and To is the pole p width. The observed resonance shape is the convolution 
of this Breit-Wigner with the photon spectrum. Figure § compares the resulting spectrum 
with the input Breit-Wigner. With the Breit-Wigner, the cross section is about 5% lower 
than for a narrow resonance with the same coupling. The phase space correction naturally 
cuts off the M nn spectrum at 2m n ; we add an upper limit M n7T < M p + 5r , matching the 



HERA analysis from which the 'jp parameters are extracted JLJ] . The correction for higher 
masses is a few percent. 

Most studies of 7^ — > pp have modelled the exclusive 7r + 7r~ spectrum with a Breit-Wigner 



p plus a non-resonant 7r + 7r term, and an interference term [18| : 



da 

IMZ, 



AJM^M P T P 



Ml, - Ml + iM p T p 



B 



(17) 



The ZEUS collaboration found A = -2.75 ± 0.04 pb 1 ' 2 and B = 1.84 ± 0.06 (pb/GeV) 



1/2 



14| . With the interference term, da/dM n7T falls rapidly at large W, and a correction for 
M n7T > M p + 5r is negligible. 

As nuclear size increases, the non-resonant term becomes less important. At lower ener- 
gies, for heavy nuclei it is at most a few percent. For large nuclei, the k dependence of B is 
unknwon, so we neglect non-resonant production. For the nuclei considered, this should be 
less than a 10% effect on the n + ir~ rate. 

The photonuclear cross sections are shown in Figure |4j. The total cross sections follow a 
similar pattern to da/dt\ t=0 . At low energies, the nuclear form factor intrudes, and eliminates 
coherent production, with the cutoff energy depending on the nuclear radius and the final 
state mass. The uj and <fi curves cross because the u includes a meson contribution that 
decreases with increasing energy, while the (ft does not and because a^N < a^N, with the 
exact relationship varying slightly with energy. 

Coherent photonuclear interactions are not a new concept. Several groups studied them 
experimentally and theoretically in the 1960s and 1970's, using real photon beams. More 
recent experiments, using virtual photons from muon scattering are less comparable because 
of the significant (compared with 1/Ra) photon virtuality. 

Most of the earlier calculations were used to extract parameters such as ayN and Ra 
from data. Our calculations are broadly similar. However, most of the earlier calculations 
used a disk model of the nucleus, with t± determined by a 2-dimensional Fourier transform 
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of a black disk, and in contributing only to a phase shift. Many of these calculations also 
included a real part for the forward scattering amplitude. For the p, the real part is 20- 
30% of the imaginary part, producing a 5% correction to the total amplitude ||. For other 
mesons, the phases should be comparable (when meson exchange is present) or smaller, 
for Pomeron-only interactions. In all cases, the correction should decrease slowly at higher 
energies. 

The experimental results can also test our calculations. Because of the complications due 
to differing (often poorly described) treatments of backgrounds, p width and the nonresonant 
7T + 7r - , we limit our comparisons to narrow resonances. Even for narrow resonances, there 
are complications. Few experiments cover the complete t range; extrapolations are required 
to find the cross sections. We compare results using da/dt\ t =o, because less extrapolation is 
required. 

An experiment at Cornell studied uo production from photons with average energy 8.2 
GeV striking a copper target, and found da/dt\ t =o = 9.6 ±1.2 mb/GeV 2 |TjJ. This compares 
well with our calculated da/dt\ t=0 = 9.5 mb/GeV 2 for that system. Another group studied 
production from 8.3 GeV photons incident on copper and lead targets, and found da/dt\ t=Q = 
4.1±0.7 and 19±3 mb/GeV 2 respectively |2(J. We predict higher cross sections, da/dt\ t=Q = 



5.4 and 38 mb/GeV 2 respectively. However, with the values of f<p and a^ used in the original 
analysis, we find da/dt\ t=0 = 3.9 and 26 mb/GeV 2 respectively, in agreement with the data. 
The discrepancy may stem from the large correction(60% for copper, 140% for lead) to cover 
the full angular acceptance. This correction is based on the optical model, and so depends 
on f^ and a^N- Still, these comparisons provide a useful check of our methods. 

III. HEAVY ION CROSS SECTIONS 

By integrating the photonuclear cross section over the photon spectrum, the total cross 
section is found: 



dN^k) _ Ui A T/ ^ _ f°° j U dN^{k) [°° da(iA -> VA) 



a(AA -> AAV) = ( dk^a^A -> VA) = f dk^ f 
J dk Jo dk Jt„ 



dt 



t=o 



(18) 
where the photon flux is from Eq. (0). The total cross sections for vector meson production 



at RHIC and LHC for a variety of beam species are shown in Table [Oil The corresponding 
production rates are shown in Table [TV]. These rates are very high, with exclusive meson 
cross sections about 10% of the total hadronic cross section for gold at RHIC, and 50% for 
lead at LHC. If the hadronic interaction model, Eq. @, is replaced with a hard sphere 
cutoff, {b > 2Ra), the rates rise about 5%. 

Previously, we considered a model where 7 A cross sections were scaled from ^p treating 
the nucleus as a black disk, with a ~ A 2 / 3 . In that model, the photonuclear cross sections 
had a different k dependence, shown by the open triangles in Fig. || This scaling predicted 
cross sections much lower than are found here |j21|. Most other studies of photon- Pomeron 



interactions have neglected the coherent Pomeron-nucleus coupling. Baur, Hencken and 
Trautman found a(Aup — > AupV)=3.5 mb using the Weizsacker-Williams photon flux for 



gold incident on a free proton target at RHIC |22 |. 



The final state rapidity is determined by the ratio of the photon energy and t\\ in the lab 
frame: 



so da/dy = kda/dk and da/dy is found by a change of variables. For the p, the cross section 
must be integrated over M nn . dN/dy is shown in Fig. |6|, for the p (Breit-Wigner), <fi and 
J/ip, for production with gold at RHIC and calcium at LHC The cross section is largest at 
low photon energies, corresponding to y < 0, because the photon flux drops as k rises. The 
peaking is larger for the p and u, where meson exchange increases the cross section at low 
energies. 

At RHIC, when photons from both nuclei are added together, separate peaks appear in 
dN/dy. At LHC, the higher energies spread the distribution over a much wider rapidity 
range, and the double peaked structure largely disappears, except for the p. 

The total perpendicular momentum spectrum is dominated by the two nuclear form 
factors, through the cutoff on the emitted photon p± and through the slope b of da/dt. The 
(quadrature) sum of these is roughly \/2?i/Ra, or about 45 MeV/c for gold. This cutoff 
is a clear signature of coherent interactions, shared with 77 and coherent double Pomeron 
interactions. Because of the exclusive production and small p±, these events can be easily 
separated from background such as grazing nuclear collisions, beam gas interactions and 



incoherent photonuclear interactions [p3 ] 



sions 



Other authors have considered coherent double-Pomeron interactions in heavy ion colli- 
p4| . If reactions where the nuclei collide are excluded, the cross section is very small, 



and depends critically on the range of the Pomeron 



IV. MULTIPLE VECTOR MESON PRODUCTION 

Because the cross sections are so large, the probability of having multiple interactions in 
a single nucleus-nucleus encounter is nonnegligible. This can be quantified by considering 
the meson production probability P(V) = da /2irbdb for a single collision. Figure [7] shows 
this probability for various mesons for (a) gold beams at RHIC and (b) lead beams at 
LHC; the probabilities reach the 1% level. For the p, we show results from two different 
hadronic interactions models: the Woods-Saxon model, Eq. (g), and also a hard sphere 
b > IRa result. Although the two models predict very similar overall rates, the probability 
is significantly affected around b = 2Ra, doubling the maximum probability at LHC. 

This probability is high enough that, even in the absence of correlations PEfl , significant 
numbers of vector meson pairs should be produced. Neglecting correlations, the pair pro- 
duction probability is P(V) 2 /2 for identical pairs and P{Vi)P(V2) for nonidentical pairs. 
These cross section and corresponding yearly rates are tabulated in Table [V]. 

Triples and higher multiples are also possible. Of course, final state interactions may 
affect what is observed in a detector. 



V. INTERPLAY BETWEEN PHOTONUCLEAR AND TWO-PHOTON 

INTERACTIONS 



Two-photon physics is expected to be a significant activity at RHIC [0 and LHC f22 



| 27| . Because the kinematics for both photonuclear interactions and two-photon collisions 
is dominated by the coherent coupling to the two nuclei, and consequently by the nuclear 
form factors, the reactions have similar kinematics. So, the specific pathway may not be 
determinable from the final state. 

The rates for photonuclear interactions are considerably higher than for comparable 
two-photon reactions. For example, the production rate with gold beams at RHIC is 
about 8 Hz. For comparison, the rates for 77 — > 7/ and 77 — > /o(980) |2l] are 0.13 and 
0.02 Hz respectively. At higher masses, the situation is similar, and photonuclear vector 
meson production may overshadow two-photon production of scalar and tensor final states, 
complicating meson spectroscopy. It is worth noting that the photonuclear pp rate is many 
times larger than that expected from 77 production |[27|| . On the other hand, the same 
techniques developed to study 77 physics at heavy ion colliders are effective at selecting 
photonuclear interactions. 

The similarity of the two interactions can lead to some interesting opportunities. For 
example, when the final states are indistinguishable, interference can occur between the pho- 
tonuclear and 77 production channels. Interference between leptonic decays of photonuclear 
produced vector mesons and Bethe-Heitler e + e~ pairs has been studied previously JT2J; at 
future heavy ion colliders, many more channels will be accessible. 

Despite the similarity, the different dN/dy spectra may allow statistical separation of 
77 and photonuclear interactions. This is because the photon flux scales as 1/k, while the 
Pomeron/meson couplings have a much weaker energy dependence. The dN/dy distribution 
for photonuclear production is broader than for two-photon production. This is especially 
true for photon-meson interactions, where the 1/k photon spectrum combines with a cross 
section that drops with energy to form two peaks at large positive and negative rapidity, 
corresponding to low energy photons from the two nuclei in turn. Figure [V| compares the 
two distributions for 77 — > / 2 (1270) and for production of a hypothetical vector meson with 
the same mass. Two vector meson models are shown: one based on the u couplings and 
another with only Pomeron coupling. The two are very similar at y = 0, but the u model 
peaks about 10% higher. 

It may also be possible to separate the two classes by their different impact parameter 
dependence. A two-photon interaction can occur at a significant distance from both nuclei, 
while a photonuclear interaction must occur inside or very near a nucleus. So, the two classes 
will have a different impact parameter dependence. In a collision, photon exchange can leave 
the nuclei excited in a giant dipole resonance, regardless of other activity; the excitation 
probability depends on the impact parameter [22|. This may allow us to statistically separate 



the two classes. 

In a broader study, it will also be possible to separate the two classes by measuring how 
the rates vary with A; the meson and Pomeron flux rises should be less A dependent than 
photons. This study will require data with a variety of different beam species. 



VI. CONCLUSIONS 

We have calculated the rates and rapidity distribution for exclusive photonuclear produc- 
tion of vector mesons in heavy ion collisions, using a Glauber model calculation. The LHC 
results involve photons with energies considerably above those currently available. Although 
this adds to the physics interest, it also introduces some uncertainty into the rates presented 
here. For the J/ip, in particular, the spectrum must soften at some photon energy, and the 
rates given here may be overestimates. 

The production rates are large enough that heavy ion colliders could be used as vector 
meson factories. The and J/ip production rates at LHC are comparable to those at existing 
or planned meson factories based on e + e~ annihilation. 

These rates are higher than comparable two-photon or double-Pomeron interactions. 
Thus, these photonuclear interactions should be the dominant source of exclusive interactions 
at heavy ion colliders. In addition to single interactions, multiple vector meson production 
from the same ion pair should be measurable. 

We thank J0rgen Randrup for suggesting the hard sphere plus Yukawa potential form 
factor. We would like to acknowledge useful conversations with Stan Brodsky and Ralph 
Engel. This work was supported by the US DOE, under contract DE-AC-03-76SF00098. 
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TABLES 



TABLE I. Luminosity and beam kinetic energy for several heavy ion beams at RHIC and LHC. 
Because RHIC will be dedicated to heavy ion acceleration, it is likely to run a wider variety of 
beams than LHC. The RHIC luminosities are from Ref. Q. Different references quote somewhat 
different heavy ion luminosity for LHC; these numbers are calculated from Table 3 of Ref. p8| , 
assuming 2 experiments, and that (as suggested) average luminosity is 45% of peak luminosity. 



Particle 



Machine 



Max Beam Energy 



Design Luminosity 



Gold 

Iodine 

Silicon 

Lead 

Calcium 



RHIC 

RHIC 

RHIC 

LHC 

LHC 



100 GeV/n 
104 GeV/n 
125 GeV/n 
2.76 TeV/n 
3.5 TeV/n 



2 x 10 26 cm- 2 s- 1 
2.7 x W^cm^s- 1 
4.4 x 10 28 ernes' 1 

1 x 10 26 cm- 2 s- 1 

2 x lO^cm- 2 ,^ 1 



TABLE II. Constants for 7p — > Vp production. The slopes by are in GeV 2 , while X and Y 
are in /xbarns, for W in GeV [11]. The f v values are from Eq. (|l0|). 



Meson 



by 



X 



Y 



ft/** 



o 



J/1> 



11 


5.0 


0.22 


26.0 


1.23 


2.02 


10 


0.55 


0.22 


18.0 


1.92 


23.1 


7 


0.34 


0.22 






13.7 


4 


0.0015 


0.80 






10.4 



TABLE III. Cross sections, in millibarns, for production of vector mesons for the beams given 
in Table 2. 



Meson 


RHIC-Au 


RHIC-I 


RHIC-Si 


LHC-Pb 


LHC-Ca 


P° 


590 


230 


8.4 


5200 


120 


LO 


59 


24 


0.9 


490 


12 


4> 


39 


14 


0.4 


460 


7.6 


JAA 


0.29 


0.11 


0.0036 


32 


0.39 



TABLE IV. Meson production rates, in Hz, at design luminosity for various beams. 



Meson 



(RHIC-Au) 



RHIC-I 



RHIC-Si 



LHC-Pb 



LHC-Ca 



o 



J/V 



120 


620 


370 


520 


230,000 


12 


63 


42 


49 


23,000 


7.9 


39 


18 


46 


15,000 


0.058 


0.30 


0.16 


3.2 


780 



12 



TABLE V. Cross sections and rates for production of vector mesons pairs in gold beams at 
RHIC and lead beams at LHC. The yearly rates are for 10 7 second/year at RHIC and 10 6 sec- 
ond/year of operation at LHC. The LHC pJ/ip rates should be used with caution because of the 
energy extrapolation. 

Meson a (fib) (RHIC-Au) Rate/year (RHIC-Au) a (fib) (LHC-Pb) Rate/year (LHC-Pb) 

LUUJ 

H> 
p°UJ 

pV 

uxfi 



720 


1,400,000 


8,800 


880,000 


6.2 


12,000 


73 


7,300 


3.8 


7,500 


76 


7,600 


133 


270,000 


1,600 


160,000 


104 


210,000 


1,640 


160,000 


9.6 


19,000 


150 


15,000 


1.3 


2,700 


200 


20,000 
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FIGURES 




1500 2000 

k [GeV] 



FIG. 1. The photon flux for gold beams using (a) the hadronic interaction probability described 
in the text, (b) the hard sphere approximation, (c), the photon flux for b > 2Ra, Eq. (^), and for 
silicon beams with the same approach as in (a). Both are at RHIC energies, with k in the target 
frame.. 
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FIG. 2. The electromagnetic form factor for gold. The solid line is the exact result from Fourier 
transformation of a Woods-Saxon potential, and the dotted line from Eq. (|14|). 
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FIG. 3. dcr/dt\t=o for coherent p, w, 4> and J/^ production on a gold target. The solid circles 
are the Glauber calculation, while the open squares and triangles show the A 2 weak absorption 
and ^4 4 ' 3 black disk scaling results respectively. The J/ip is in the weak interaction regime, while 
the lighter mesons are closer to the black disk. The Glauber calculation removes much of the 7p 
cross section energy dependence visible with the simple scaling. 
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FIG. 4. Photonuclear cross sections a(jA — ► VA) for p° (solid line), ui (dashed line), 
line) and J/ip (dot-dashed line) coherent production on a gold target. 
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FIG. 5. Predicted M^ spectrum for exclusive p photoproduction (solid line), compared with 
the input Breit-Wigner spectrum (dashed line). 
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FIG. 6. The rapidity distribution of produced p (top), <fi (center) and J ftp (bottom). The left 
panels show are for gold beams at RHIC, while the right panels are for calcium beams at LHC. 
The solid line is the total production, while the dashed line is for photons coming from the nucleus 
on the left. 
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FIG. 7. Probability of vector meson production as a function of collision impact parameter for 
(a) gold at RHIC and (b) lead at LHC. The solid lines show the Woods-Saxon hadronic interaction 
probability, while the dashed lines shows the result of a hard sphere calculation for the p. 
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FIG. 8. Comparison of the rapidity distributions for two-photon production of the /2(1270) 
(dashed histogram) and two models of photonuclear production of a hypothetical X(1270), for 
gold beams at RHIC. The solid histogram is for an X(1270) with couplings matching the w, while 
the dotted histogram corresponds to Pomeron coupling only. 
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